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Introduction 

Interactions among the atmosphere, terrestrial ecosystems and 

yearTflthough m h oT fh ‘ he S “ bJeC ‘ ° f ,nvesti S atl °" «t many 
(Hewlett and Hibberl 1967: WoSTa ^gs^The t^'lT'h' 
including the effects of climate and hydro.og^ on veXtarionThe 

cvrl S °[j egetati ° n ° n h y droI °gy- the effects of the hydrological 
ycle on the atmosphere, and interactions of the cycles Ada material 
Dux such as solutes and trace gases. The intent of this pape^Ts to 
denhly areas of critical uncertainty, discuss modeling approaches 

lng I^msfd “ Pr ° blems ' and then P r °Pose techniques for test- 
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" septate 60 “ ““‘'‘r 1 ”''' ^ (3) terrestrtal carbon cycle, 
but for r* . lss ues of process and scale somewhat artificially 
but for convenience In discussing the Issues more clearly Issues „f 

undeMo"d "h " 6 b ‘° l0glCal or P h y steal Processes are not well 

msponle to mil?! ’ ‘ h ' b ‘° loglcal C ° ntr0,s over ^system 
esponse to C0 2 fertilization are not known, although hypotheses 

abound. Issues of scale are familiar in both the physica/Lences 

where constderaWe problems with cloud parameteriLtlons persTs!’ 

and biology, where the problem of extrapolating from the organism 

level (where our understanding ,s eoncentndedf to regmns 3 ^ 

globe ,s unresolved. In some cases, the Issues ofproeL anTscaS 
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cNaughton s (1986) classic paper on scaling from leaf to canopy: 
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“rencrterre^ria'lecolo^ at the global scale, yet C0 2 exchange 
is mediated by single leaves and within subsys- 
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m derfaWHty torepUcate system behavior along an environmental 
models ability to repiica y anolication of these tech- 

*7on^e sm.uiat.ons (Hobbs et ah. In press). Brether. 
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^ , i tVi^t thev embody our best understanding. 

SirS St- ^ experiment, by theoty or 
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to their interest and challenge. 
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production and the outcome of plant competition. Several factors 
are involved in this plant-radiation interaction. First, photosynthesis 
(A) is highly and nonlinearly sensitive to incoming photosyntheti- 
cally active radiation (PAR), with a response curve often fitted by a 
hyperbola (Johnson and Thomley, 1984). 

Changes between cloud and clear sky radiation will have sub- 
stantial effects on photosynthetic rates and associated evapotran- 
spiration (ET). Because these rates Eire highly nonlinear, the para- 
meterization chosen for clouds in a coupled atmosphere-biosphere 
model will have large consequences during climate change simula- 
tions. That is, taking a grid cell with humidity greater than thresh- 
old for cloud formation and assigning an average cloudiness will 
result in an average radiation field. Because the photosynthesis-PAR 
curve is strongly nonlinear, carbon and water exchange for that grid 
cell will be calculated incorrectly, possibly resulting in subsequent 
errors in the atmospheric water. 

At steady state, this problem may be resolved by implicitly para- 
meterizing the model to take into account “typical” cloud statistics, 
and this is often done. For nonsteady-state simulations or those of 
trace gas-modified climate, a physically based cloud statistic for 
each grid box is more desirable. The response of the vegetation can 
be correctly modeled by integrating the A-PAR function over the 
cloud distribution field for each grid cell. 

The above argument assumes that the A-PAR relationship adjusts 
instantaneously and that PAR is the sole control over A. In fact, 
plant responses include lags and are controlled by multiple factors 
(Knapp and Smith, 1988; Schimel et al., 1991). First, while the drop 
in A following a decrease in PAR is rapid, the increase in A after PAR 
is increased may be slow (Knapp and Smith, 1990). Thus, an 
increase in cloudiness could produce a larger decrease in A and 
evapotranspiration than predicted with a linear extrapolation as a 
result of the hysteresis of the A-PAR relationship in time. This effect 
may be small at the global scale but has not been evaluated. 

In water-stressed environments, cloudiness may actually 
increase net carbon gain. Water-stressed vegetation at high PAR 
will exhibit near-total stomatal closure and high resistance (g) to 
water and C0 2 exchange. Thus in dry environments and under 
clear sky conditions, A may be near zero and ET low for much of 
the day, despite high photosynthetic capacity at full sunlight when 
well watered. Cloud passage under such conditions may allow 
stomatal opening and permit some photosynthesis. Knapp and 
Smith (1988, 1990) have shown that typical levels of cloudiness 
increase net daily A over that predicted for clear sky conditions in a 
subalpine environment. 
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This opens the following questions: (1) To what extent are such 
responses (hysteresis of A-PAR) tuned to the environment by evolu- 
tionary adaptation, or (2) how and how fast can plants respond to a 
changing light intensity? (3) Is the response physiological, in which 
case we may assume that a near-optimal response will keep pace 
with the changing environment, or (4) is the response genotypic, in 
which case the response will occur over time scales of community 
and evolutionary change (Field, 1990; Schimel et al., 1990). 

In addition to the physiological processes, light also affects the 
outcome of competition between plant species. The widely used 
Jabowa family of models (e.g., Pastor and Post, 1986) includes the 
effects of light competition as a primary factor structuring forest 
communities. Each species included in the model has parameters 
describing its tolerance of high and low light levels, and shading is 
computed using a simple radiation model. While the response of 
ecosystem models to changing temperature and water has been 
evaluated (Schimel et al„ 1990, 1991; Clark, this volume; Pastor 
and Post, 1988), effects of changing light environment have not been 
simulated. Indeed, scenarios of changing earth surface radiation 
from GCMs have not been widely released and may be quite unreli- 
able (Smith and Vonder Haar, 1991). In fact, plant responses to the 
absolute light environment (as opposed to the relative light environ- 
ment defined by shading) may not be well known. 

Returning to the issue of cloud parameterizations in atmospheric 
GCMs (AGCMs), it is clear that simulated cloudiness structure or 
statistics in the AGCM not only must satisfy the requirements of 
correct simulation of the cloud radiation feedback on mean global 
and regional temperature. It also must accurately model the vari- 
ability of clouds and hence the solar radiation field at the earth’s 
surface for correct simulation of plant response. The time scales of 
these two requirements may well be different. For example, the 
earth radiation budget is normally calculated with a time step of 
about a month, and radiation calculations In most AGCMs are inte- 
grated at longer time steps than many other processes. For simula- 
tion of A, ET, and g, statistics describing the variability in monthly 
cloudiness but for morning, midday, and afternoon conditions may 
be required. Cloud parameterizations in AGCMs should be devel- 
oped with requirements based on sensitivity analyses of ecosystem 
models as one criterion. In addition, the use of cloud statistics in 
ecosystem simulations should be investigated further. 

' One type of model subsystem checking is conceptual checking to 
insure that no important process has been omitted. The above dis- 
cussion points out a series of interactions that have not been con- 
sidered in the development of most models of land surface interac- 
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tions. Assuming that the above interactions were to be included in 
an earth system model, how could they be tested? This question has 
a number of answers. First, the physiological response of plants to 
cloudiness can be investigated using leaf and canopy measures of 
gas exchange, and some such studies are under way (Knapp, per- 
sonal communication; Knapp and Smith, 1988, 1990). In addition, 
many data sets collected using eddy correlation measurements of 
water and C0 2 exchange include incoming radiation or even cloud 
statistics (e.g., the first ISLSCP [International Satellite Land Surface 
Climatology Program] Field Experiment [FIFE] Information System, 
NASA-Goddard Space Flight Center) and could be analyzed for 
effects of cloudiness and hysteresis. Cloud passage is often treated 
as a noise term but could be analyzed as a signal. 

The physiological adaptation of plants to new cloud regimes could 
be simulated experimentally using growth facilities or transplants 
but will be difficult to relate to the field situation and to global 
change. The effects of changing light environment on plant competi- 
tion can be modeled using improvements of extant models but will 
be well-nigh impossible to test in the field given the time scale of 
plant succession. In this area, sensitivity analysis of rigorously 
developed models may have to suffice, possibly augmented with field 
studies using analogous gradients. 


Soil Moisture 


Soil moisture is a key parameter linking atmosphere and bios- 
phere. Soil moisture is a key control over decomposition (Parton et 
al., 1987), over leaching of nutrients, and over plant growth. The 
central role of soil moisture storage and of runoff to the biota and 
hydrology make their correct representation in earth system models 
crucial. In order to simulate the surface energy balance, climate 
models use the basic prognostic equation for soil moisture (w): 


9(w) 

at 


r - ET + m - f 


( 1 ) 


where r is rainfall, ET is evapotranspiration, m is snowmelt, and f is 
runoff. Recent work has focused on improving the representation of 
ET (Sellers et al., 1986). This approach is acceptable for simulation 
of sensible and latent heat exchange over broad areas, but is it ade- 
quate for linking climate, plant growth, and biogeochemistiy? 

There are several issues to consider. First, plant available water 
(pw) is defined by the following equation in real landscapes: 

8(pw) r 

^ — f t,om It — + m - f + i -d] 


( 2 ) 
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where pw is plant available water; t is time; om is organic matter 
content; r, ET, m, and f are as above; d is deep drainage; and i is 
run-on. The function must be converted from total water (the con- 
served quantity) to plant available water using the soil moisture 
release function, which varies largely as a function of soil texture 
and soil organic matter content. In many ecosystems, run-on is a 
significant control over landscape-scale plant production as much of 
the production may be concentrated in zones of run-on (Schimel et 
al., 1985; Noy-Meir, 1977). Similarly, trace gas efflux may be much 
greater from zones of water concentration than from upland areas 
(Schimel et al., 1985; Parton et al., 1988), both because of higher 
microbial activity and because of the erosional concentration of 
nutrients in run-on zones (Schimel et al., 1986). 

The above issues raise problems of scale. The simple prognostic 
equation for soil water, as modified to include detailed representa- 
tion of evapotranspiration, currently suffices to simulate atmos- 
pheric moisture and surface energy balance in AGCMs. These repre- 
sentations may also suffice for calculations of water balance for 
large drainages (Gleick, 1987). When the focus of interest changes 
to the biosphere, these representations become less satisfactory. 
Stored soil water is an influence over plant production in many 
ecosystems, and this is of course influenced by position in the 
runoff/run-on continuum. Field studies and simulations of nitrous 
oxide (N 2 0) and ammonia (NH 3 ) flux show position in the landscape 
along this continuum to be of predominant importance (Schimel et 
al 1986; Parton et al., 1988). Similarly, methane flux is influenced 
by meso- and microscale hydrology and topography (Whalen and 
Reeburgh, 1988). The above problems require mapped or statistical 
data on topography, soil properties, vegetation, land use, and engi- 
neering structures for resolution in models. Clearly, some of these 
problems require resolution higher than achievable or desirable in 
an earth system model; equally clearly, these processes (productiv- 
ity, trace gas flux) are of significance to the earth system. 

Improved coupling of atmospheric and hydrological models is also 
important to understanding the interactive role of fresh water 
ecosystems in the earth system. Changes in climate will have signifi- 
cant consequences for rivers and lakes, and these changes will 
undoubtedly have feedbacks, at least to regional climate and cer- 
tainly to human societies. Current AGCMs include no representa- 
tion of fresh water interactions and barely even include them due to 
coarse resolution of geography. Yet fresh water systems are critical 
to nutrient and organic matter transport, interact strongly with 
marginal wetlands, and are the dominant vector for nutrient trans- 
port to the oceans. There have been repeated critiques of the lack of 


‘David J>. SchimeC 


411 


consideration of coastal and shelf regions of the oceans in global 
carbon cycle models; simulation of these areas will require consider- 
ation of the inputs from fresh water systems. Indeed, as Broecker 
(1981) has pointed out, fresh water may play a dominant role in mil- 
lennial changes in ocean circulation and climate. 

A critical link in the land surface-atmosphere-hydrology part of 
the “wiring diagram” then is the interface between global and 
regional changes in P-E, and their expression in real hydrological 
systems. A variety of approaches has been used to develop such 
interfaces. Wilks (1988) used a statistical conversion from AGCM- 
predicted changes to “real” weather and soil properties from a crude 
data base to translate between AGCMs and detailed plant models 
but did not address run-on /runoff and other landscape processes. 
Schimel et al. (1990) used a similar approach for biogeochemical 
cycling. Gleick (1987) used AGCM output as input to a basin hydrol- 
ogy model, but focused on water output rather than on distribution 
and effects within a watershed. Giorgi (1990) used a mesoscale 
atmospheric model forced by AGCM output and a land surface para- 
meterization to examine effects of mountainous topography and 
sub-AGCM grid circulations; this study did not consider routing of 
runoff. Vordsmarty et al. (1989) used a basin-scale hydrological 
model to study the routing of runoff within a large drainage (the 
Amazon) and its effects on evapotranspiration and nutrient flux. 
This model could be forced by simulations such as those of Dickin- 
son (1987) or Shukla et al. (1990) to evaluate regional hydrological 
consequences of a changed vegetation in the Amazon region. 

The problems of including hydrology in earth system models are 
serious, however. Wood et al. (1990) have pointed out that most 
hydrological models are extensively parameterized to account for the 
peculiar features of topography, soils, and geology in a catchment 
area of interest. This is because of the complex nonlinear interac- 
tions that occur within catchments, and the detailed geographic 
data required to analyze those interactions. The response of catch- 
ments is also very sensitive to initial conditions (Wood et al., 1990; 
Figure 1), again resulting in complex responses. Work by Wood, 
Band, Gupta, and others (Wood et al., 1990; Wood et al., 1988; 
Gupta et al., 1986; Band and Wood, 1988) has attempted to develop 
scale-independent measures that can be used to develop simplified 
models of watershed response, at least above threshold sizes [Repre- 
sentative Elementary Areas: Wood et al., 1990). These theoretical 
approaches and simulations hold promise for developing improved 
representation of land surface hydrology and its interactions with 
biology and the atmosphere in the near future. Improvement of the 
representation of hydrology in earth system models will greatly aid 
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Figure 1 . Results of simulations of runoff on a Kansas tallgrass prairie water- 
shed (Konza Prairie Research Natural Area Long-Term Ecological Research Site; 
also, the site of FIFE. Adapted from Wood et al, 1990). Patterns are superim- 
posed on the digital elevation model of the site, (a) Rainfall for a storm on August 
4, 1 987; (b) predicted runoff assuming dry initial conditions. 

in improving the biogeochemistry and biology in these models, but 
improvements to the runoff routing will be required in addition to 
improved ET simulations. 

Validation of earth system models including hydrology is again 
difficult. Replication is essentially impossible at the scales of inter- 
est and experimental manipulation unlikely. Two approaches have 
been used with some success. The large held experiment, exempli- 
fied by FIFE or the Hydrological Atmospheric Pilot Experiment (Sell- 
ers et al., 1988], can test parameterizations of the ET and small- 
scale routing issues by extensive measurements of soil moisture, 
stream flow, and fluxes above the canopy. Basin-scale investigations 
can compare model simulations to river hydrographs (Vorosmarty et 
al., 1989), although the validity of the internal processes remains 
untested. In the end, these models will have to be tested against a 
few intensive field studies and long-term records and then compared 
broadly to routinely collected data on rainfall, stream and river flow, 
and river biogeochemistry. Simulations of paleoevents, such as the 
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Figure 1(b). 


hypothesised Younger Diyas event, will also be helpful in establish- 
ing confidence in these complex simulations. Finally, the global 
parameterizations should be based on extensively tested and more 
complex regional models. 

The Carbon Cycle 

Issues in the terrestrial carbon cycle have been extensively 
argued over the past decade or so, and many unresolved questions 
remain (Houghton et al., 1983; Detwiler and Hall, 1988; Tans et al., 
1990). Tremendous effort has been expended on resolving the ques- 
tion of whether the terrestrial biosphere is a source or sink for CO 2 . 
I will not attempt a new analysis of this question but will instead 
outline steps that seem necessary to resolving it, and further steps 
required for inclusion of an interactive terrestrial biosphere in earth 
system models. The questions that must be answered to address the 
source/sink issue seem now to bee (1) Where and when, if ever, is 
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C0 2 fertilization expressed? (2) What are the geography and carbon 
budgets for disturbed ecosystems; can these systems be a sink? (3) 
Can climate change (variability) over the past century have caused a 
directional change in carbon storage? In order to include an interac- 
tive biosphere into earth system models, the following additional 
questions must be answered: (1) If terrestrial ecosystems are now a 
sink, what are the limits to sink strength? (2) How will vegetation 
distribution change and how will this affect the carbon cycle? 
Finally, for both issues (source/sink, predictive models), how can 
the answers be tested at large scales? 

Can the terrestrial biosphere serve as a long-term sink for car- 
bon? Enhanced C0 2 will lead to increased photosynthesis under 
certain conditions. Three factors may limit the significance of this 
increase. First, it may be a transient, as some evidence shows that 
plants may acclimate to increased C0 2 and show a gradually 
decreasing response. 

Second, the effects of C0 2 may be attenuated by constraints from 
other limiting factors. Nutrients or water may limit C0 2 uptake at lev- 
els only slightly greater than under current C0 2 concentrations. The 
interaction of C0 2 fertilization with other limiting factors requires far 
more study in a range of ecosystems with varying limiting factors. 
C0 2 fertilization has the potential to alter nitrogen and water use by 
allowing increased enzyme efficiencies (photosynthetic engines con- 
tain large amounts of N) or by increasing water-use efficiency through 
stomatal effects. Under most circumstances, these changes in water 
or N-use efficiencies will result in the production of plant tissue with 
reduced content of N and other nutrients. This is because C0 2 fertil 
ization does not enhance nutrient availability through any known 
mechanism and, with increased efficiency, more biomass will be pro- 
duced on the same amount of nutrients. The production of plant tis- 
sue with higher C: element ratios will increase microbial uptake of 
nutrient when that tissue is decomposed, competing with plants and 
reducing nutrient availability (Parton et al., 1987; Melillo et al., 1984). 
This feedback will tend to reduce the effects of C0 2 fertilization on pri- 
mary production, homeostatically. We have evaluated the conse- 
quences of increasing C0 2 on carbon storage in a grassland model 
(Figure 2) and found the behavior described above to apply; that is, 
feedback through plant-microbial competition for nutrients limited 
the effects of C0 2 -induced increases in resource use efficiencies. 

Third, an extrapolation of the above effect suggests that the abil- 
ity of terrestrial ecosystems to store carbon is limited by the avail- 
ability of other resources. If C0 2 were to Increase carbon storage in 
some or most ecosystems, how much carbon could they store? A 
certain amount of C can be stored in plant biomass, especially in 
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wood. However, the maximum biomass of forests is constrained 
mechanically, by nutrients and light, and by human requirements 
for forest products. Carbon storage in soil organic matter is a more 
permanent sink, having a turnover time of hundreds to thousands 
of years. All organic matter, however, has nutrients associated with 
it, and soil organic matter is normally quite rich in nutrients. For 
example, N in humus is generally 5 to 20 times more concentrated 
than in plant biomass (Paxton et al., 1987). Thus any major increase 
in soil carbon will constitute a major sink for N and other nutrients, 
eventually leading to the same sorts of restrictions on primary pro- 
ductivity described above. In effect, the humus becomes a competi- 
tive sink for nutrients with the plants, and production must decline. 
Bogs and other wetland areas where peat accumulates are the con- 
temporary examples of this process, having high carbon storage and 
low primary productivity. 

The use of fertilization in agriculture to support or enhance the 
direct C0 2 effects on primary production would seem to enhance a 
negative feedback to atmospheric C0 2 concentration. However, rates 
of decomposition are relatively high in agricultural systems and 
rates of carbon stabilization low. Most agricultural soils contain at 
the most a few percent carbon, and the carbon content tends to 
decline slowly or stabilize at low levels (Bouwman, 1990). Thus the 
negative feedback due to fertilizing agricultural soils is limited; how- 
ever, carbon once stabilized in cultivated soils may have a long resi- 
dence time (Parton et al., 1987; Trumbore et al., 1990). 

Disturbed ecosystems may pass through a period of rapid carbon 
accumulation, where stocks of biomass and soil carbon are replen- 
ished (see Houghton et al., 1983). This phenomenon is evident in 
recently harvested forests, after burning, and following severe storm 
damage. The expansion of the use of terrestrial ecosystems in the 
past century may have “reset the clock” for many ecosystems to the 
point where they are accumulating organic matter rapidly and serv- 
ing as a sink. While the accumulation of carbon following distur- 
bance is included explicitly in the budgetary estimates of Houghton 
et al. (1983) and Detwiler and Hall (1988), our data on the geogra- 
phy of disturbance and recovery rates are not numerous (Botkin 
and Simpson, 1990). 

In conclusion, terrestrial ecosystems are unlikely to be a long- 
term sink for C0 2 although certain areas act as sinks for some 
intervals. The ability of ecosystems to store carbon is limited and 
unlikely to accommodate the extent of atmospheric C0 2 increases 
projected as a consequence of industrial civilization. 

If greenhouse gas emissions lead to increased temperatures, this 
will also affect the carbon cycle. Increased temperatures, other fac- 
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tors being equal, will accelerate decomposition and cause loss of 
stored carbon, increasing the atmospheric inputs of C0 2 (Schimel et 
al., 1990). In general, decomposition is quite sensitive to tempera- 
ture (e.g., Parton et al., 1987). In a recent simulation exercise, we 
showed that across a range of grasslands, increased temperature 
resulted in soil carbon losses despite enhanced production due to 
C0 2 enrichment (Figure 2). While this result is only for one ecosys- 
tem type, other models exhibit similar behavior for other ecosystem 
types (Pastor and Post, 1986), suggesting it may be a fairly general 
result. If this Is true, then C0 2 fertilization responses may not keep 
pace with respiration losses resulting from increases in tempera- 
ture. It is also possible that climate changes over the past few cen- 
turies, including a possible greenhouse effect term, have resulted in 
significant changes In soil carbon storage. Substantial improve- 
ments in data and models will be required to test this hypothesis 
and even to establish the direction of the change; nonetheless, this 
would be an interesting exercise even with current process models. 

Testing process-based models of the carbon cycle presents major 
challenges. While regional models can be developed and tested effec- 
tively, repeating this process in every ecosystem globally is pre- 
cluded by logistics. Testing at the global scale is feasible using 
atmospheric C0 2 concentration gradients, but these tests do not 
test the internal dynamics of the model very well, and hence the 
predictive capability of the model remains in question. 

I suggest a strategy based on three techniques. First, models 
should be based on understanding of principles that have been 
proven to be broadly applicable — a theory-based approach — rather 
than on empirical relationships (Schimel et al., in press). Real 
progress is being made in this area for both autotrophic and het- 
erotrophic processes, although empirical models remain the domi- 
nant type for global representations. 

Second, while traditional field tests cannot be made everywhere, 
remote sensing can obtain data in virtually all terrestrial ecosys- 
tems. Remote observations can serve to quantify temporal and spa- 
tial distributions of ecosystem properties with current techniques 
(Tucker et al., 1986). Excellent progress is being made on retrieving 
more quantitative information on hydrology, physiology, and biogeo- 
chemistry using satellite data (Schimel et al.. 1991b; Wessman et 
al., 1988). While many of the traditional formulations for ecological 
and hydrological processes are not readily tested with remote sens- 
ing, often these models may be transformed or reparameterized to 
allow such testing (Schimel et al., 1991b). 

Finally, the aggregate output of models already tested regionally 
using remote sensing and field techniques can be compared to 
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